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ABSTRACT 

We model the broad-band X-ray spectrum of Cyg X-3 in all states displayed by this source 
as observed by the Rossi X-ray Timing Explorer. From our models, we derive for the first 
time unabsorbed spectral shapes and luminosities for the full range of spectral states. We 
interpret the unabsorbed spectra in terms of Comptonization by a hybrid electron distribution 
and strong Compton reflection. We study the spectral evolution and compare with other black 
hole as well as neutron star sources. We show that a neutron star accretor is not consistent with 
the spectral evolution as a function of Le and especially not with the transition to a hard state. 
Our results point to the compact object in Cyg X-3 being a massive, ~ 30Mq black hole. 

Key words: gamma rays: observations - radiation mechanisms: non-thermal - stars: individ- 
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1 INTRODUCTION 

Cygnus X-3 is one of the brightest as well as one of the first known 
X-ray binary systems. It was discovered already 40 years ago (Gi- 
acconi et al. 1967) and has been observed by all X-ray missions 
as well as extensively observed in other wavelengths from radio to 
PeV energies, although without any confirmed detection above the 
hard X-ray regime. An interpretation of its intrinsic X-ray proper- 
ties and spectral evolution is, however, still missing from the liter- 
ature and the source remains enigmatic. Despite several recent at- 
tempts to apply physical models to explain its intrinsic broadband 
spectrum (Vilhu et al. 2003; Hjalmarsdotter et al. 2004; Hjalmars- 
dotter et al. 2008) including a detailed modelling of the wind pa- 
rameters (Szostek & Zdziarski 2008), the underlying spectral vari- 
ability and the true accretion geometry have yet to be explained. 

In this paper, we make use of recent constraints from Hjal- 
marsdotter et al. (2008), hereafter Hj08 and Szostek & Zdziarski 
(2008), hereafter SZ08, to interpret the intrinsic spectral variability 
observed in Cyg X-3 and discuss its spectral evolution in compar- 
ison with other X-ray binaries. For the first time, we present unab- 
sorbed spectral shapes and luminosities characteristic for the full 
spectral variability of this source. We use archival data from the 
Rossi X-ray Timing Explorer (RXTE) which include the full range 
of observed spectral states. 
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Figure 1. The absorbed spectral shapes of Cyg X-3 from Szostek & 
Zdziarski (2004). The 42 observations are averaged into 5 groups, the hard 
state (blue solid line), the intermediate state (cyan long dashes), the very 
high state (magenta short dashes), the soft non-thermal state (green dot- 
dashes) and the ultrasoft state (red dotted line). 
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2 THE RXTE DATA 

We use the data from 42 pointed observations with the RXTE Pro- 
portional Counter Array (PCA) and the High Energy X-ray Tran- 
sient Experiment (HEXTE) between 1996-2000 for which PCA 
data from Proportional Counter Units 0-2 were available. The PCA 
spectra were extracted in the 3-25 keV range using only the top 
layer of the detectors. The corresponding HEXTE spectra were ex- 
tracted in the 15-110 keV energy band, using both HEXTE clus- 
ters. We allow a free relative normalization of the HEXTE spectra 
with respect to that of the PCA. The observations are listed in Ta- 
ble 1. The data were previously published in Szostek & Zdziarski 
(2004). They found that the 42 observations could be divided into 
five groups based on the shape of the observed absorbed spectrum. 
Unabsorbed spectral shapes and luminosities were however not dis- 
cussed in that preliminary study and it was not clear whether the 
observed differences were due to absorption effects or to intrinsic 
spectral evolution. The data is also the same as used for the clas- 
sification of X-ray and radio states in Szostek, Zdziarski & Mc- 
Collough (2008). All the observations and their division into the 
five groups are shown here in Fig. Q] The figure replaces fig. 1 in 
Szostek & Zdziarski (2004) and in Szostek et al. (2008) where there 
is a mismatch between colours and linestyles. 

Here, we re-model the average spectra of each group for a de- 
tailed investigation of the physics of the observed spectral states, 
as well as the 42 individual spectra to construct colour-colour and 
colour-luminosity diagrams. A 1 per cent systematic error was 
added to each of the averaged spectra. We label the spectra of group 
1-5 in Szostek & Zdziarski (2004) and Szostek et al. (2008), in or- 
der of increasing softness, as the hard state, the intermediate state, 
the very high state, the soft non-thermal state and the ultrasoft state, 
respectively. 



3 SPECTRAL MODEL 
3.1 Continuum model 

We use a Comptonization model, eqpair (Coppi 1992, 1999), in 
which soft seed photons are up-scattered in a hot flow located inside 
the inner radius of or on top of an accretion disc. The source of the 
soft seed photons is assumed to be an accretion disc with the spec- 
trum of a disc-blackbody (diskbb, Mitsuda et al. 1984, included 
in eqpair). The slope of the Comptonized spectrum is determined 
by the ratio of heating of the Comptonizing electrons (apart from 
Coloumb interactions and Compton heating) to cooling due to the 
injection of soft seed photons. This can be expressed as the ratio 
of the hard and soft compactnesses 4/4 where £ = Lo"T/(rm e c 3 ), 
L is the luminosity, <r T is the Thomson cross section, r is the char- 
acteristic size of the X-ray emitting region and m c is the electron 
mass. The model is weakly dependent on the value of f s , through 
Coulomb interaction and pair production only, and thus 4 is poorly 
constrained in the fits. We therefore freeze it at a value of 4 = 100, 
corresponding to a high luminosity and a small radius of the Comp- 
tonizing plasma, appropriate for the case of Cyg X-3, and similar 
to the case of GRS 1915+105 (Hannikainen et al. 2005; Zdziarski 
et al. 2005). The plasma optical depth, t 1o1 , includes a contribution 
from electrons formed by ionization of the atoms in the plasma, r e , 
as well as a contribution from e ± pairs, r tot - r e . The electron distri- 
bution, including the temperature, T c , is calculated self-consistently 
from energy balance and can be purely thermal or hybrid if an ac- 
celeration process is present. The acceleration rate is assumed to 
have a power-law shape with index T for Lorentz factors between 



Table 1. The log of the RXTE observations. The spectral classes 1, 2, 3, 4 
and 5 correspond to the observations (22-23, 25-30), (1-8, 24), (9-13, 19, 
21), (41^12), and (14-18, 20, 31-40), respectively. 
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1.3 and 1000. In the case of a hybrid plasma 4 = 4h + 4th, i-e., the 
thermal heating plus the power supplied to accelerated electrons. 
This model for the continuum intrinsic spectrum is the same as used 
by Vilhu et al. 2003, Hjalmarsdotter et al. 2004, Hj08 and SZ08, 
and its application to Cyg X-3 is discussed more elaborately espe- 
cially in Hj08. The model is also described in detail by Gierlinski 
etal. (1999). 

Compton reflection is included (pexriv, Magdziarz & 
Zdziarski 1995), parametrized by a solid angle R = Q/2n. In an un- 
obscured geometry R corresponds to the true solid angle of the re- 
flector as seen from the primary X-ray source. If the primary source 
is obscured, R also involves the relative amount of emission seen in 
reflection to that observed directly, and can thus greatly exceed 2, 
the maximum value for an isotropic source in an un-obscured situ- 
ation (see further discussion in Hj08 and SZ08). The strength of the 
reflected component also depends on the system inclination, which 
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is not well known in Cyg X-3. The strong orbital modulation sug- 
gests a high inclination, but not too high since there is no eclipse 
by the secondary. For the purpose of modelling the reflection, we 
use 60° (Vilhu et al. 2007). We allow the reflecting medium to be 
ionized but find that our fits are consistent with a neutral reflecting 
medium. 

3.2 Absorption effects in the wind 

Absorption in Cyg X-3 is a complex matter. Studies with e.g. Chan- 
dra (Paerels et al. 2000) have revealed a rich line complex with 
both absorption and emission features at energies <10 keV. SZ08 
recently studied the effect of the stellar wind from the Wolf-Rayet 
companion on the spectra of Cyg X-3 using medium resolution 
BeppoSAX data and modelled the wind using the photoionization 
code cloudy (Ferland et al. 1998). They found that the emission 
and absorption structure at low energies required a two-phase wind 
made up by a hot tenuous plasma together with cool dense clumps 
filling ~ 1 per cent of the wind volume. Since the RXTE data 
lack high spectral resolution and cover only energies above 3 keV, 
and since we are more interested here in the intrinsic spectra than 
the wind parameters, we use a simplified model for the absorption 
based on the findings of SZ08. 

In a first attempt aiming to reproduce the physical conditions 
of the wind as well as possible, although in a simplified approach, 
we model the wind using absori (Done et al. 1992), with a weakly 
ionized one-temperature plasma with initial values for the optical 
depth, temperature and ionization parameter set to average values 
from the results in SZ08. We find that freezing the values for the 
temperature and ionization parameter to the average values from 
SZ08 does not give satisfactory fits to any of the data sets, show- 
ing that a one-phase single-temperature and single-ionization ab- 
sorber is not a good approximation to the complex wind structure. 
However, allowing column density, temperature and ionization free 
in our fits, we find that our data can not constrain these values or 
break the degeneracy between them. In most cases, the obtained 
error ranges are consistent with neutral absorption. 

On the other hand, we find that, just like in Vilhu et al. 2003, 
Hjalmarsdotter et al. 2004 and Hj08, very good fits to the data can 
be obtained using a two-fold neutral absorber (absnd) with one 
medium with column density A'He.tot, covering the entire source and 
one medium with column density A'He.par, covering a fraction f c , of 
the source (with the column densities expressed in He rather than 
H, given that Wolf-Rayet stars contain very little hydrogen, see fur- 
ther discussion on abundances below). This model was also used by 
SZ08 as a first approximation for the effect of the wind absorption 
in order to calculate the ionizing spectrum for the wind modelling. 
The overall similarity of the obtained intrinsic spectral shapes in 
that first step to those obtained in the final wind modelling in SZ08, 
shows that this two-fold model, although more phenomenological 
than physical, gives a very good approximation to the effect of the 
wind absorption in Cyg X-3. In this paper, where we aim at finding 
physically consistent models for the intrinsic spectral shapes rather 
than any detailed modelling of the absorbing medium, we therefore 
use the neutral two-fold absorber to model the continuum absorp- 
tion since it gives a better phenomenological approximation of the 
true wind effects than any simplified ionized absorber. We find that 
for our purposes, this is a sufficiently accurate approach, but note 
that the values of the column density can not really be constrained 
based on our data. 

We first fit the data with the abundances of a Wolf-Rayet star 
of WNE type from Hamann & Grafener (2004), the same abun- 



dances as used by SZ08. Since our simplified modelling of the 
absorber does not take into account ionization, it is, however, not 
consistent with high iron abundances for column densities close to 
the ones fitted by SZ08. Furthermore, we find that for some spec- 
tral states where reflection plays an important role, especially in 
the hard state (see Section 4.1) the RXTE data do not seem to be 
consistent with high iron and metal abundances. We thus use nor- 
mal solar abundances (Anders & Ebihara 1982), but with the metal, 
A >Hi; (only in eqpair), and iron, A Fe , abundances (requiring the 
same value for the absorption as in eqpair) as free parameters. 
The abundances are thus allowed to differ between different spec- 
tral states. This does not correspond to a physical reality, but is 
rather a compensation of not treating ionization, which may well 
differ between the different spectral states. We note that an anoma- 
lously low iron abundance, much lower than the values in our best 
fit models, were in fact suggested by Terasawa & Nakamura (1994, 
1995). We give column densities of the local absorber as re-scaled 
to helium as Af He . 

In Hj08, an additional Thomson absorber was used to account 
for additional attenuation caused by electron scattering in the stel- 
lar wind, presumably responsible for the strong orbital modulation. 
The phase averaged data was corrected by increasing the normal- 
ization by 1/0.75 times to roughly correspond to phase maximum. 
Here, we reconsider this assumption. If the wind is highly ionized 
scattering does not remove photons so the photons missing at phase 
minimum have to show up somewhere else. The simplest hypothe- 
sis is that they get back-scattered, and show up around phase max- 
imum. This means that when we look at the source at the phase 
where the optical depth is lowest, we see photons back-scattered 
from the direction where the optical depth is the highest. Averag- 
ing over the phase, we see all the photons, and no correction is 
required. We do not know the exact degree of ionization of the 
scattering medium. The wind-modelling in SZ08 suggests, how- 
ever, that the hot tenous highly ionized phase of the wind makes 
up 99 % of its volume. We thus assume that photons are not sig- 
nificantly removed in the scattering process and do not make any 
corrections to the derived luminosity due to scattering. This gives 
the most conservative estimate for the luminosity and we note that 
the actual luminosity may be higher. 

Finally, we fix the interstellar absorber to Mi.gai = 1.5 x 10 22 
cm" 2 (e.g. Chu & Bieging 1973). 

3.3 The iron complex 

A strong edge at 7 keV apparent in both data sets modelled in 
SZ08 implies the presence of neutral iron, not accounted for in their 
model of the wind. This confirms the presence of at least some ad- 
ditional cool material in the system, either in the form of a reflector, 
or as an absorber in addition to the Wolf-Rayet wind, or both. In our 
models to the RXTE data, the edge results from a combination of 
neutral absorption and reflection. An additional edge at ~9 keV is 
also required by the data. In the wind modelling by SZ08, this edge 
arises from highly ionized iron in the otherwise optically thin hot 
inter-clump wind medium. For simplicity, we model this hot phase 
with an edge with energy £ e( ] ge and optical depth T edge . We find our 
values of the optical depth of this edge to correspond well to the 
column density of the hot phase in SZ08. 

The presence of a range of ionization species of iron in the Cyg 
X-3 wind also gives rise to a rich emission line spectrum, resolved 
by Chandra (Paerels et al. 2000) and partly resolved by BeppoSAX 
(SZ08). The iron line complex is also clearly visible in the RXTE 
data, but due to the low spectral resolution of the PCA instrument, 
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a single broad gaussian centered at ~ 6.5 keV is sufficient to model 
it. For a more physical approach, however, we separate the most 
prominent wind emission line at 6.7 keV, as observed by Chan- 
dra and BeppoSAX, from any possible additional flourescent line 
arising from either absorption or reflection. In the case of strong 
reflection with R > 1 an accompanying flourescent line is indeed 
expected. We thus model the iron line complex with two gaussians, 
one with a fixed energy E = 6.7 keV and fixed physical width 
cr = 0.1 (corresponding to an arbitrary thin line with the PC A res- 
olution) and allow an additional line with energy E Fc between 6.4 
and 7.5 keV, physical width o"f c < 0.4 keV and equivalence width 
EWpe We find that including even more line structure, e.g. a line at 
6.9 keV, also seen in higher resolution data, is not required in our 
fits to the RXTE data. 

We note that the large error ranges on most of our best-fit pa- 
rameters in Table 2 indicate that they are not at all very well con- 
strained. Some parameters, e.g. the column density Nu c , can simply 
not be constrained by our data. For several of the other parameters, 
however, the large error range is due to the presence of several, 
themselves much more narrow, local minima. Including all these 
local minima results in a very large global error range. We have 
found that freezing the blackbody temperature at the best-fit tem- 
perature when determining the error ranges for the other parameters 
reduces the number of local minima, especially those with intrin- 
sic spectral shapes and luminosities far from the best fit model, and 
especially in the hard and intermediate states. This procedure was 
followed for all states except the ultrasoft state. The problem with 
very large error ranges and unconstrained column densities, how- 
ever, remains. We also can not fully rule out the presence of addi- 
tional minima in this very complicated multidimensional parameter 
space. 



4 SPECTRAL STATES 
4.1 The hard state 

The hardest state of Cyg X-3 differs from the canonical hard state, 
as seen in e.g. Cyg X-l, in that it peaks at much lower energies, 
~ 20 keV, compared to other black holes as well as neutron stars 
in their hard states. Such a low cut-off implies an electron tem- 
perature of the Comptonizing electrons below 10 keV, usually ob- 
served only in soft states where strong cooling is provided by the 
massive influx of soft seed photons from a strong blackbody com- 
ponent. In addition, it was found by Hj08 that any attempt to fit 
the data with thermal Comptonization proved difficult, since the 
high-energy slope does not match well a cut-off due to the finite 
temperature of a Comptonizing electron distribution (similar to but 
not equal to an e-folded power-law). Hj08 showed that the Cyg X-3 
hard state could indeed be fitted with a model with a high ampli- 
fication factor (4/4 > 7), despite the low cut-off energy, if the 
power to the electrons in the plasma was provided entirely by ac- 
celeration. A model where the apparent hard state was instead only 
a result of increased absorption in the stellar wind also gave a good 
fit to the INTEGRAL data in Hj08. The best fit in that study was 
found for a model dominated by Compton reflection, where most 
of the direct emission was obscured from view. 

The RXTE data of Cyg X-3 in its hard state used here consist 
of the spectra from 8 observations, made between July 1997 and 
August 1998. The variations between them are small and the indi- 
vidual spectra as well as the averaged one is similar to the hard state 
as observed by INTEGRAL in Hjalmarsdotter et al. (2004) and in 



Hj08. For the sake of statistics, we derive our best fit models based 
on the averaged spectrum of all 8 observations. 

We find that all three models tested in Hj08 give good fits 
(X 2 /y < 1-2) also to the RXTE data, with the reflection dominated 
model again providing the best fit (y 2 /v = 0.48) (We note that our 
very low x 2 values for the RXTE data may well be due to an over- 
estimation of the systematic errors added to our data). Despite its 
ability to fit the data, this model was however ruled out as a phys- 
ical description of the hard state of Cyg X-3 by Hj08, since it, just 
like the wind absorption model, presented in the same paper, re- 
quired a very luminous soft incident spectrum. This is not in agree- 
ment with the other results of that study which strongly preferred 
a 'real' state transition rather than obscuration effects (see further 
discussion in Section 6.1). In SZ08, BeppoSAX data from a hard 
state was modelled with a continuum model with equally strong 
reflection (R = 10) as in the reflection model of Hj08 but with a 
harder incident spectrum with a very high blackbody temperature 
and almost completely non-thermal distribution of the Comptoniz- 
ing electrons. Such a spectrum does not give a satisfactory fit to the 
RXTE data wifh;f 2 /v ~ 2. 

However, we find a very good fit (x 2 /v = 0.44) for a model 
with moderately strong reflection and with an incident spectrum of 
moderate hardness and luminosity. The model is plotted together 
with the data in Fig. [2^, and the components of the unabsorbed 
model are shown in Fig. [2p. The best fit parameters are listed in 
Table 2. 

In our model, the observed spectrum above 10 keV is domi- 
nated by strong Compton reflection that creates the characteristic 
shape of the cut-off around 20 keV. The model requires R = 2.5, 
which corresponds to a situation just exceeding maximum possible 
reflection in a geometry where most of the direct emission is still 
observed directly (for i = 60°). 

The incident spectrum is that of a moderately hard state. With 
4/4=2.0, the power supplied directly to the electrons is higher 
than the power in the irradiating seed photons, but due to the low 
electron temperature of only 12 keV, the Comptonized spectrum 
(the sum of the thermal and non-thermal Comptonized compo- 
nents) before reflection is flat in vFv space. The best-fit blackbody 
temperature of 0.6 keV suggests that the disc is truncated at some 
radius far from the innermost stable orbit, but not as far away as 
in the canonical hard state of e.g. Cyg X-l (Gierlinski et al. 1997), 
or as suggested by the non-thermal model in Hj08. We note, how- 
ever, that the exact blackbody temperature, in our proposed model 
located at energies much lower than the RXTE bandpass at 3 keV, 
is not well constrained by the data. 

In addition to the Fe xxv line at 6.7 arising in the wind, an ad- 
ditional broader line centered at 6.4 keV is required by the data, in 
agreement with strong reflection from a neutral medium. We note 
that the fitted EW^ of the broad line is almost twice as high as 
expected (calculated self-consistently using thcompfe Zdziarski, 
Johnson & Magdziarz 1996; Zycki, Done & Smith 1999) from the 
fitted R using the best fit abundances. Using the WR-abundances 
from Hamann & Grafener (2004), however, the fitted line EWp e 
agrees well with the expected EWp e for both the hard and the inter- 
mediate state. 

It is clear that even if the accretion geometry is similar to that 
of the canonical hard state, the physics of this accretion flow dif- 
fers from that of a standard optically thin advective flow where t ~ 1 
and the electron temperature is observed to be between 50-100 keV 
for black holes (e.g. Zdziarski & Gierlinski 2004), and 40-60 keV 
for neutron stars (Barret 2001). The optical depth of r = 4.8 in- 
cludes a small fraction of pairs. A 38 per cent non-thermal fraction 
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Figure 2. a: Averaged data from eight observations with the RXTE/PCA (magenta) and RXTE/HEXTE (blue) of Cyg X-3 in its hardest state together with 
the model described in Section 4.1, absorbed (dashed) and corrected for absorption (solid line), b: The components of the unabsorbed model: the unscattered 
blackbody (magenta long dashes), Compton scattering from thermal (green short dashes) and non-thermal (blue dot-dashes) electrons, reflection and the Fe K 
line (red dotted line) and the Fe 6.7 keV emission line from the wind (cyan dotted line). 




Figure 3. a: Averaged data from 9 observations by RXTE/PCA (magenta) and RXTE/HEXTE (blue) of Cyg X-3 in its intermediate state, together with the 
model described in Section 4.2, absorbed (dashes) and corrected for absorption (solid line), b: The components of the unabsorbed model: the unscattered 
blackbody (magenta long dashes), Compton scattering from thermal (green short dashes) and non-thermal (blue dot-dashes) electrons, reflection and the Fe K 
line (red dotted line) and the Fe emission line at 6.7 keV (cyan dotted line) from the wind. 



of the electron distribution with injection with power-law index 3.9 
is required to fit the highest energy HEXTE data. We note that the 
HEXTE spectrum above 100 keV does not agree with that observed 
in the hard state by SPI in Hj08 (fig. 8 in that paper) but agrees 
with the ISGRI data (same figure). According to the model pre- 
sented here, Cyg X-3 is not likely to be detected by GLAST in its 



hard state, unless the maximum Lorentz factor of the accelerated 
electrons is larger than the assumed 1000. 

We note that a slightly softer version of our moderate reflec- 
tion model gives a rather good fit to the BeppoSAX data of SZ08, 
and a slightly better fit than the phenomenological model presented 
there. This, however, only holds if the abundances are allowed to, 
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as here, differ from normal Wolf-Rayet ones, which they were not 
in SZ08, and which may not, as discussed in Section 3.1, corre- 
spond to a physical reality. We have also found this model to give 
a good fit to the INTEGRAL data, better than the wind-absorption 
and the non-thermal models in Hj08, but not as good as the reflec- 
tion model presented there. The INTEGRAL data are indeed much 
more difficult to model in part due to the large discrepancy between 
the SPI and ISGRI data in Hj08. Since the model presented here for 
the hard state is also agreement with the results from Hj08 in that it 
represents a state with the Compton amplification factor 4/4 > 1, 
and probable truncation of the inner disc, we suggest it as a better 
physical description of the intrinsic spectrum of Cyg X-3 in its hard 
state than any earlier presented models. 

Using our best-fit model of the averaged data as an input, 
we also model the eight individual spectra to derive their intrinsic 
shapes and colours (see Section 5). Due to uncertainty of whether 
the small variations between the individual spectra is due to intrin- 
sic spectral variability or to a change in the absorbing medium, we, 
however, refrain from discussing and/or interpreting this intra-state 
variability in terms of change in the parameter values. The same 
procedure is performed for all five spectral states. 



4.2 The intermediate state 

An intermediate state generally refers to a spectral state in between 
the hard and the soft state and is usually observed in transitions. 
In Cyg X-3, state transitions are frequent and an intermediate state 
should therefore be rather common. The RXTE data presented here 
consist of 9 spectra, one group of continuous observations in Au- 
gust 1996, judging from the ASM lightcurve just before a transition 
from a hard to a soft state, plus one observation from September 
1996 after a short transition into a hard state and back to a soft. 

Our best fit model to the Cyg X-3 intermediate state is similar 
to that of the hard state, with a similar blackbody temperature, but 
with heating of the electrons balanced by cooling by seed photons 
to give 4 /4s = 0.96. The main difference in the observed spectral 
shape is the somewhat lower amplitude of reflection than in the 
hard state. A harder injection spectrum of the non-thermal electrons 
also gives more power at energies above a few hundred keV. The 
bolometric X-ray luminosity is somewhat higher than in the hard 
state. The model is shown together with the data in Fig. [3^, and its 
components in Fig. [3p. The best fit parameter values are listed in 
Table 2. 

We note that our model predicts a strong annihilation feature. 
The actual presence of such a feature (at 5 1 1 keV, above the fitted 
range here) depends on the unknown compactness (t s ), and will be 
weaker at lower compactnesses. In addition, a lower value of the 
maximum Lorentz factor of the accelerated electrons would pro- 
duce less high-energy photons and thus less pairs at a given com- 
pactness, thus reducing the feature. 

4.3 The very high state 

The term 'very high state' is used in the literature to describe a va- 
riety of spectral shapes, not without some confusion. It usually de- 
scribes a spectral state with the presence of both a strong blackbody 
component and a strong component from thermal Comptonization, 
typically with the Comptonized component starting already at the 
peak of the blackbody (e.g. Done & Gierlihski 2004). In some 
sources the spectral shape in this state is often more or less identical 
to the intermediate state (as pointed out by Zdziarski & Gierlihski 



2004), and differs only in amplitude. The state may or may not be 
the most luminous one. 

The very high state as displayed by Cyg X-3 is a rather soft 
variant of this state with a strong Comptonized component, rather 
different from its intermediate state. This very high state is, judg- 
ing from both the RXTE/ ASM lightcurve (high flux and moder- 
ate hardness) and pointed observations by both RXTE and INTE- 
GRAL, the second most common spectral state in Cyg X-3 after 
the hard state (Hj08). Fig. [4^ shows our best-fit model together 
with the data averaged from 7 observations by RXTE in Febru- 
ary, June and July 1997. The spectrum is similar to that in Vilhu 
et al. (2003) but modelled here in a somewhat different way. In 
our interpretation, its unabsorbed spectral shape is similar to that 
of the softer states of GRS 1915+105 and with similar parameters 
(e.g. Zdziarski et al. 2005). The blackbody temperature is 1.0 keV, 
and a strong disc component is present in the spectrum. Cooling 
is strong with 4/4 = 0.33 but the spectrum is still dominated by 
the Comptonized component. The plasma has a temperature of 5 
keV and an optical depth of 3.4. Reflection is strong also in this 
state, R = 2.3, but since the spectrum is much softer, reflection 
does not have as strong an effect on the observed spectral shape. As 
opposed to the hard and intermediate state, in this very high state 
the broad iron line has an EWp c somewhat lower than that expected 
from fluorescence corresponding to the strength of the continuum 
reflection component (for both the fitted abundances and for typical 
WR-abundances). The RXTE data of Cyg X-3 in its very high state 
is consistent with a rather steep distribution of non-thermal elec- 
trons and less emission in the model at energies > 100 keV. The 
very high state is the most luminous state in Cyg X-3, calculated 
from the absorbed as well as the absorption corrected bolometric 
X-ray flux. The inferred geometry is a disc extending close to the 
innermost stable orbit with the hot Comptonizing flow formed on 
top of the disc rather than within its inner radius. 



4.4 The soft non-thermal state 

Cyg X-3 also displays a soft state with a moderate disc component 
and a very strong nonthermal tail. This state was observed in April 
2000, right after a major radio flare. The spectral state is similar to 
that labelled 'high' state in Done & Gierlihski (2004, fig. la), and 
appears to be a transitional state between the ultrasoft and the very 
high states. In terms of luminosity, however, this state in Cyg X-3 
is equal to that of the intermediate state and thus the least luminous 
of the observed soft states in this source. 

Our modelling of this state includes a Comptonized blackbody 
similar to that of the very high state, but weaker, plus a signifi- 
cant non-thermal component requiring an almost 75 per cent non- 
thermal distribution of the Comptonizing electrons. The model is 
shown in Fig. [5] The best fitting parameters are given in Table 2. 



4.5 The ultrasoft state 

The softest state of Cyg X-3 belongs to the so called ultrasoft state, 
a state with the soft component dominating the spectrum and with 
only weak emission above ~20 keV. Such a spectrum is believed to 
arise at high accretion rates where the thin disc extends all the way 
into the innermost stable orbit of the black hole (or neutron star). A 
strong non-thermal tail, with no apparent cut-off below 200 keV, is 
also present. The ultrasoft state is generally observed at accretion 
rates ranging from about 0. 1L E up to L E or higher (e.g. Gierlihski & 
Done 2004) while being absent in sources not reaching luminosities 
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Figure 4. a: Averaged data from 7 observations by RXTE/PCA (magenta) and RXTEfHEXTE (blue) of Cyg X-3 in its very high state, together with the model 
described in Section 4.3, absorbed (dashes) and corrected for absorption (solid line), b: The components of the unabsorbed model: the unscattered blackbody 
(magenta long dashes), Compton scattering from thermal (green short dashes) and non-thermal (blue dot-dashes) electrons, reflection and the Fe K line (red 
dotted line), and the Fe 6.7 keV emission line from the wind (cyan dotted line). 




Figure 5. a; Averaged data from 2 observations by RXTE/PCA (magenta) and RXTE/HEXTE (blue) of Cyg X-3 in the state labelled the soft non-thermal state 
together with the model described in Section 4.4, absorbed (dashes) and corrected for absorption (solid line), b: The components of the unabsorbed model: the 
unscattered blackbody (magenta long dashes), Compton scattering from thermal (green short dashes) and non-thermal (blue short dashes) electrons, reflection 
and the Fe K line (red dotted line) and the Fe 6.7 emission line from the wind (cyan dotted line). 



close to 0.1L E , like e.g. Cyg X-l. The state is in some sources the 
most luminous one. 

RXTE has observed Cyg X-3 in its ultrasoft state on several oc- 
casions, typically right before a major radio flare. Here we model 
the spectrum averaged from a total of 16 observations in June 1997, 
February 2000 and June 2000. Fig. [6p, shows the data together 



with the best fit model with parameters in Table 2. In our unab- 
sorbed version, the spectrum of Cyg X-3 is indeed very similar to 
the ultrasoft spectra of e.g. GRS 1915+105, XTE J1550-564 and 
GX 339-4 (see a comparison of the absorbed spectral states in fig. 
8 in Zdziarski & Gierlinski 2004). Fig. [6p shows the components 
of the model. The spectral de-composition shows that the Cyg X-3 
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Figure 6. a: Averaged data from 16 observations by RXTE/PCA (magenta) and RXTE/HEXTE (blue) of Cyg X-3 in its ultrasoft state, together with the model 
described in Section 4.5, absorbed (dashes) and corrected for absorption (solid line), b: The components of the unabsorbed model: the unscattered blackbody 
(magenta long dashes), Compton scattering from thermal (green short dashes) and non-thermal (blue dot-dashes) electrons, the Fe 6.7 line (cyan dotted line) 
and an additional broad Fe line (red dotted line). 



spectrum in this state is not really blackbody dominated, but dom- 
inated by the component arising from thermal Comptonization of 
the strong blackbody by electrons of moderate temperature. This 
is similar to the ultrasoft state of GRS 1915+105 in Zdziarski et 
al. (2005). Just like in GRS 1915+105, simple modelling with a 
blackbody + power-law model of this state would thus lead to an 
overestimation of the blackbody temperature. Our best-fit model 
includes a nearly flat non-thermal tail. If such a tail is real it may 
be strong enough to be detected by GLASTfLKT after one year of 
operations (according to the GLASTfLAH performance web pag^3)- 

Our model is consistent with zero reflection. Models with high 
reflection, R ~ 5, are, however, found within A% 2 just above 3.0. 
We note that the BeppoSAX data modelled in SZ08 do not seem 
to be consistent with zero reflection, at least not without a neutral 
absorbing medium in addition to the ionized wind in their model. 
The highly reflection dominated model presented in SZ08 is, how- 
ever, problematic since the strong reflection is not accompanied by 
any fluorescent Fe line apparent neither in their BeppoSAX data, 
nor seen in e.g. Chandra data in the Cyg X-3 ultrasoft state (M. 
McCollough, in prep.) In the model presented here, it seems on the 
other hand as if a broadened line, in addition to the narrow line at 
6.7 keV is required by the data, despite the zero reflection, unless 
the 6.7 keV line itself is allowed to be broadened. We conclude that 
due to the low resolution of the PCA instrument, the data used here 
do not allow for any detailed line diagnostics, but that our zero- 
reflection model for the continuum in the ultrasoft state agrees with 
the Chandra data. 

The unabsorbed spectral shapes according to our best fit mod- 
els of the full range of the observed spectral variability are shown 
in Fig.|7]together with the envisaged acccretion geometries. 



www-glast.slac.stanford.edu/software/IS/glast_lat_performance.htm 



5 SPECTRAL EVOLUTION 

5.1 The colour-colour diagram of Cyg X-3 

Having derived unabsorbed spectral shapes of the range of spec- 
tral states displayed by Cyg X-3, we can now calculate its intrinsic 
X-ray 'colours' and study the spectral evolution in a colour-colour 
diagram. To enable comparison with other classes of X-ray bina- 
ries, black holes as well as neutron stars, we use the same colours 
as Done & Gierlinski (2003, hereafter DG03), where the soft colour 
is defined as the ratio of the intrinsic fluxes in the 4—6.4 keV and 
3-4 keV bands, and the hard colour as the ratio of the intrinsic 
fluxes in the 9.7-16 keV and 6.4-9.7 keV bands. In Fig. [8ji the 
colour-colour diagram of Cyg X-3 is plotted together with the data 
of the black holes and atoll sources from DG03, Gladstone, Done 

6 Gierlinski (2007) and Done, Gierlinski & Kubota (2007). Note 
that the colours of Cyg X-3 are not corrected for reflection, i.e. they 
represent the observed although absorption corrected spectral shape 
and are therefore not dependent on the actual spectral decomposi- 
tion of the unabsorbed spectra. 

It is obvious from Fig.[8]that the track in the colour-colour di- 
agram traced out by Cyg X-3 does not resemble the characteristic 
Z-track of the atoll sources (the Z-like shape is actually more pro- 
nounced in the atolls than in the Z-sources, DG03; no Z-sources 
are included in Fig. [8). It is rather more similar to that of the black 
holes, except that the the high values of the soft colour of Cyg X-3 
in its soft states keep it from entering the lower left corner of the 
diagram, the area occupied only by black holes and inaccessible 
to neutron stars as defined by DG03. The high values of the soft 
colour in Cyg X-3 is due to the high blackbody temperature that re- 
sults in a hard slope of the spectrum below 5 keV in the soft states. 
The soft states of Cyg X-3 occupy similar colours to those of GRS 
1915+105, but it is obvious that Cyg X-3 spans a much wider range 
of hard colours than does GRS 1915+105, which does not show 
transitions into a hard state due to its constantly high luminosity, al- 
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Table 2. The parameters of the best-fit models, as presented in Sections 4.1^1.5. The uncertainties correspond to 90 per cent confidence, i.e. A^ 2 = +2.71. The 
large error bars are a result of the presence of several themselves much more narrow minima. 



hard state intermediate state very high state soft non-thermal ultrasoft state 



iV He ,i 
fi 

A Fe 

A>He 
^edge 
^"edge 

kT* 

4/4 
4th/ 4 
r b . 

rnj 
Te 
C 

R 



il/2n 



EW Fe 

fbol,abs 
F g 

bol 

Hoi 

x 2 /y 



10 22 cm 
10 22 cm 



keV 



keV 



keV 
keV 
eV 

erg cm~ 2 s -1 
erg cm -2 s 
erg s 



4 +0 - 4 
1.7+ 1 - 7 


7 +al 
u -'-0.2 

6 +al 
q 37+O.I6 

-0.03 

0.6+ 03 
2 oo+8'06 

dsM 

O+0.1 
-0.1 
5 +0.12 
°-0.14 

4.77 
8.3 

6.4+ - 2 
3 +0 ' 

"'■'-0.3 

310 
7.7 x 10-' 
1.6 x 10~ 8 
1.6 x 10 38 
47/93 



0.17+ 



3.' 
4.76+ 



4 +a2 
o.^_ 04 

UJ -0.4 

9 +0 3 

7 +a2 
u '-0.1 

9 28if£ 



0.25 
0.6; 
0.90 



12 



n to+0.49 

O c+0.6 

T 7<: + 1.01 
j.;j„ 39 

3.78 



7.6 

22+25 

64 ;8: 5 2 

0.4_ .2 
190 
7.9 x 10 
2.3 x 10 
2.2 x 10 
23/84 



6 +0 2 
1.9™ 



0.5+ 



S3 

-0.3 

0.8+°? 
1.4$ 
9,19+0.10 

21^ 

-0.03 

1.0+° 4 

36 + ^ 9 

-0.05 

36 +! > W 

—0.19 

T C+0A 

J - -0.7 
t 40+0.14 

3.40 
4.6 

7 O + 1.0 

Z - J -1.7 

6 5 +0 - 2 
0.4_ .2 
120 
1.1 x 10~ 8 
4.9 x 10~ 8 
4.7 x 10 38 
24/87 



6 +0 - 2 

, 9 +3.5 

4 +o: ^ 

u -0.2 

4 +a2 
u -^-0.2 

1.0/ 

9 14+0.09 

71 -0.10 

0.27 



1.1 



70-05 



0.32 
0.71 



m 

-0.09 
9 c+0.4 
Z - J -0.1 

4 70+ 4 - 72 
4 -' u -1.51 

4.74 
2.2 



2.0 + 



0+5.9 

6.5+ - 1 

2+^ 

-0.2 

200 
8.2 x 10~ 9 
2.5 x Kr 8 
2.5 x 10 38 
45/74 



0.< 



7 +o.i 

0.2+ - 1 
6+^ 

3+0.6 
-0.5 
1.0/ 

9.14+0.07 
34+^ :o ' 

u -0.04 

1.2+0-! 
0.21 
0.50+! 



01 



B 

0.02 

2 l +al 
' 4,1 -0.1 
4 72+0-68 
' -0.26 

5.15 
2.8 



0.0 



+0.1 

6 5+ 01 
"■ J -o.i 

0.4_ .i 
270 
8.7 x 10" 9 
3.7 X Kr 8 
3.6 x 10 38 
60/89 
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b The acceleration power-law index T for the non-thermal electrons was limited between 2.0 and 4.0. 

c Calculated from the energy balance, i.e. not a free fit parameter. 

d The energy of the additional broad iron line was limited between 6.4-7.5 keV. 

c The physical width of the additional broad iron lime was limited to < 0.4. 

8 The bolometric flux of the unabsorbed model spectrum normalized to the PCA data. 

h The unabsorbed bolometric luminosity assuming a distance of 9 kpc (Dickey 1983, assuming 8 kpc for the distance to the Galactic Centre; Predehl et al. 
2000). 



ways exceeding 0.3L E (Done, Wardziriski & Gierlinski 2004, here- 
after DWG04). Overall, the track of Cyg X-3 in the colour-colour 
diagram is rather similar to that of GRO J1655-40 (DG03, fig. 2), a 
black hole transient harbouring a 7M© black hole (Orosz & Bailyn 
1997) and, just like GRS 1915+105 and Cyg X-3, a strong radio 
source. 



5.2 The colour-luminosity diagram of Cyg X-3 

The mass of the compact object in Cyg X-3 is unknown. To study 
the spectral evolution as a function of L E , we thus have to assume 
a probable mass for the compact object. Alternatively, we may use 
the track traced out by Cyg X-3 in its colour-colour and colour- 
luminosity diagram to by comparison to other objects make an ed- 
ucated guess about the approximate mass and nature of its X-ray 
emitting object. 

The Eddington luminosity for the X-ray emitting compact ob- 
ject is 

L E = 4nfi e GMm p c/(TT, (1) 

where fi e = 2/(1 + X) is the mean electron molecular weight, X 
is the hydrogen mass fraction, and cr T is the Thomson cross sec- 
tion. Assuming, despite the simplifications made in modelling of 
the absorption, that the accreted matter contains no or very little 



hydrogen, L E is that of a helium-star, L E = 2.5 x 10 38 M/M© erg 
s" 1 . 

In Fig. [8p the hard colour is plotted against the luminosity 
from Table 2 as L/L E for three different assumptions for the mass 
of the compact object in Cyg X-3. Only the results from the mod- 
elling of the average spectral shapes are plotted for clarity. Over- 
plotted is the hard colour versus luminosity in Eddington units for 
the black holes and neutron stars (atolls) from DG03, Gladstone 
et al. (2007) and Done et al. (2007) (the same data as in Fig. [Hi, 
for comparison. (Note that the lines connecting the points for the 
different mass-assumptions do not suggest the source evolution but 
are simply there to guide the eye). 

The track to the right assumes M\ = 1.4Mq, representing a 
neutron star accretor. If the compact object is a neutron star, its path 
traced out in the CL-diagram differs significantly from that of the 
normal atoll sources. The luminosity is close to and sometimes ex- 
ceeds L E . This behaviour is not observed in atolls but seems to be 
common in several Z-sources, neutron star X-ray binaries with con- 
stant high accretion rates (see e.g. DG03). As an extreme example, 
the Z-source Cir X- 1 can reach highly superEddington luminosities 
and spans a range of < 10L E (DG03). These sources that are con- 
stantly accreting at high accretion rates do not, however, show state 
transitions. Their tracks in a colour-colour diagram span only a nar- 
row range in hard colours, corresponding to a variable black-body 
temperature, and never displays hard colours exceeding 1.0, corre- 
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Figure 7. a: Spectral shapes of Cyg X-3 according to our models, corrected for absorption. The hard state (blue solid line), the intermediate state (cyan long 
dashes), the very high state (magenta short dashes), the soft non-thermal state (green dot-dashes) and the ultrasoft state (red dotted line), b: Proposed accretion 
geometries including radio behaviour. The non-thermal electron contribution is shown as blue stars. 



sponding to the realm of hard states in the colour-colour diagram 
(see e.g. DG03, fig. 8). Thus, the fact that Cyg X-3 shows state 
transitions is incompatible with its compact object being a neutron 
star since state transitions to the hard state is not compatible with a 
constant high accretion rate. 

The middle track is that for a compact object of M = IOMq, 
a typical mass for Galactic black holes. With this assumption, the 
luminosity span of Cyg X-3 corresponds to -0.06-0. 2L E , compara- 
ble to some of the more luminous black hole systems, like eg. XTE 
J1550-564 and LMC X-3 (DG03, fig. 2). The transition to the hard 
state would however have to take place at ~ 0. 1L E . Such high tran- 
sition luminosities are indeed observed in transient systems (e.g. 
0.3L E for GX 339-4, Zdziarski et al. 2004), despite the fact that 
theoretical models of the advective solutions that are believed to 
describe the accretion flow in the hard state are stable only up to 
luminosities of a few per cent L E (e.g. Esin et al. 1997). The fact 
that these high transition luminosities are observed only in tran- 
sient systems, and only in the hard-to-soft transition, point to them 
being a result of the transient behaviour. When a transient enters 
an outburst, it starts from quiescence. If the local accretion rate and 
thereby the luminosity increases faster than there is time to build up 
an accretion disc (usually believed to take place on the viscous time 
scale), the source may stay in a hard state up to a higher luminosity 
than what would have been possible if an accretion disc was al- 
ways present and the transition from a dominant hot advective to a 
cool disc solution was set solely by the properties of the respective 
solutions, irrespective of previous behaviour. The transient sources 
thus show hysteresis, i.e. the hard-to-soft transitions take place at 
higher luminosities than the soft-to-hard transitions. Cyg X-3, just 
like Cyg X-l, does not show any hysteresis. The transitions be- 
tween the hard to soft and the soft to hard state always take place at 
the same luminosity. 

For the transition to take place at a luminosity corresponding 



to 0.03 L E which is the case in e.g. Cyg X-l and LMC X-3 (DG03), 
both persistent systems like Cyg X-3, the mass of the compact ob- 
ject in Cyg X-3 would have to be as high as 30Mq (for L to equal 
0.03L E in its intermediate state). The track in the colour-luminosity 
diagram for Cyg X-3 assuming a compact object of M = 3OM0 is 
shown to the left in Fig. [Hp. With this assumption, the luminosity 
span in Cyg X-3 corresponds to 0.02-0. 07L E , comparable to that 
of Cyg X-l. The full colour-luminosity diagram for M = 30Mq, 
including all 42 observations is shown in Fig. [9] The pattern is 
somewhat similar to that of LMC X-l and GRO J1655 (DG03). 
Note that part of the spread in luminosity in the individual obser- 
vations is due to orbital modulation. The difference in luminosity 
between observations at phase maximum and minimum differs by 
approximately a factor of 2. 



We would like to point out that the values for the luminosity in 
Fig. [Hp are all lower limits. For the states where reflection exceeds 
2.0, the implied intrinsic luminosity should be multiplied by a fac- 
tor of minimum R/2 (assuming full sky coverage of the reflector 
as seen by the emitter) and maximum R (for a slab-geometry) to 
account for the part of the emission obscured from view. (See fur- 
ther discussion in Hj08 and SZ08). If scattering in the stellar wind 
results in a net decrease of photons (as assumed in Hj08 and SZ08 
but not here, see Section 3.2), this would also increase the intrinsic 
luminosity. The higher luminosity, the higher mass is needed to ex- 
plain the data and our value for the required mass of the compact 
object of 30Mq is thus the most conservative one. 
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Figure 8. a: Colour-colour diagram of Cyg X-3. All 42 individual observations are shown as black triangles overplotted on the data from the Galactic black 
holes (cyan circles) and atolls (yellow filled circles) from DG03, Gladstone et al. (2007) and Done et al. (2007), and GRS 1915+105 (red circles) from 
DWG04. b: Colour-luminosity diagram of Cyg X-3 assuming a mass of its compact object of 1.4, 10 and 30 Mq, from right to left (black triangles, only 
average points are shown for clarity). Data from the same Galactic black hole sources (cyan circles), atolls (yellow filled circles) and GRS 1915+105 (red 
circles) are overplotted. 




Figure 9. The full colour-luminosity diagram for Cyg X-3 assuming a black 
hole mass of 30Mq overplotted on the same data from other black hole and 
atoll-sources as in Fig 8. About a factor of two in the spread in luminosity 
is due to orbital modulation. 



6 DISCUSSION 

6.1 State transitions and the presence of a hard state in Cyg 
X-3 

An important question regarding the spectral variability in Cyg X-3 
has been whether it actually displays a state transition from a soft 
disc dominated state to a hard state with disc truncation in the same 
way as most other black hole binaries, or whether the observed 
spectral variability at energies < 10 keV is caused by variable local 
absorption. The absorbed spectral shapes resemble those of other 
black hole binaries (Szostek & Zdziarski 2004), with the exception 
of the hard state that peaks at considerably lower energies than what 
is usually observed, as discussed in Section 4.1. The similarity of 
the spectral shape of the hard state at energies >20 keV, together 
with the presence of strong absorption, affecting all states up to 
higher energies than usual for normal Galactic absorption columns, 
may at first seem suggestive that the hard state is just a strongly 
absorbed version of one of the softer states (e.g. Hjalmarsdotter et 
al. 2004). Understanding the spectral variability in Cyg X-3 is thus 
very much dependent on understanding the nature of the hard state 
and assessing whether the transition to it marks a 'real' state tran- 
sition or not. The issue of the nature of the hard state of Cyg X-3 
was investigated thoroughly by Hj08. Results from a study of the 
X-ray lightcurve showed a bimodal distribution of the soft X-rays. 
This, together with the anticorrellation of soft and hard X-rays as 
well as between soft X-ray flux and hardness, is suggestive of a 
bimodal behaviour of the accretion flow, with intrinsic spectral piv- 
oting rather than variable absorption being responsible for the ob- 
served state transitions. The similarity of the radio/X-ray correla- 
tion pattern in the hard state of Cyg X-3 to that of other sources 
in the hard state further suggested that the accretion flow in the 
Cyg X-3 hard state shares some common physics and/or geometry 
inherent to the canonical hard state. The fact that the orbital modu- 
lation was found to be stronger in the soft state than in the hard also 
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speaks against the hard state being an artefact of increased wind 
absorption. These results led Hj08 to conclude that the hard state 
of Cyg X-3 is likely to be a lower accretion rate state with a trun- 
cated or absent inner disc and that a state transition indeed takes 
place in Cyg X-3 and represents an intrinsic re-configuration of the 
accretion flow. 



6.2 The location of the inner disc 

Together with the uncertainties about the exact properties of the 
absorber, the true blackbody temperature (and thereby the inner 
disc radius), especially in the hard state is hard to determine by 
spectral modelling. We would like to stress that this is not just a 
question of the high bandpass (>3 keV) of the RXTE data. The sit- 
uation is almost equally problematic when using data extending to 
lower energies. When modelling the BeppoSAX data in SZ08, the 
blackbody temperature in the hard state could only be determined 
to within A^ 2 /v of 0.14 for a range of 0.4-2.8 keV with the same 
phenomenological treatment of absorption as used here, and had to 
be frozen when modelling the wind parameters. The high black- 
body temperature found as the best fit in SZ08 (1.7 keV) may well, 
as discussed in that paper, correspond not to the inner disc temper- 
ature but to some soft excess, often present in the spectra of black 
hole binaries (eg Cyg X-l, Gierlinski et al. 1997). In this paper, 
we find a model for the hard state consistent with a lower black- 
body temperature (0.6 keV).The exact black body temperature in 
the hard state remains, however, poorly known. 

6.3 Compton reflection and the geometry of the accretion 
flow 

A reflection dominated spectrum in the hard state of Cyg X-3 was 
first suggested by Hj08 as the best fit model to the INTEGRAL data 
presented there. Also SZ08 found that a reflected dominated spec- 
trum gave the best fit to BeppoSAX data of Cyg X-3 in both its 
hard state and a state similar to the ultrasoft state presented here. 
In this paper, our best fitting models to four of the five observed 
states (all except the ultrasoft state) require strong reflection, ex- 
ceeding R = 2.0 (for ( = 60°), the maximum solid angle in an 
unobscured geometry. This implies that if the inclination is > 60° 
the direct emission from the inner accretion flow is to some extent 
obscured from view, even if the models presented here do not imply 
luminosity corrections for the obscured parts larger than the true lu- 
minosity being 1.25 times the observed one (as compared to a cor- 
rection factor >5 in SZ08). Reflection domination is commonly ob- 
served in obscured AGN (e.g., Matt, Guainazzi & Maiolino 2003). 
Recently, a Compton reflection dominated spectrum has also been 
suggested to explain a low non pulsating state of the neutron star 
binary GX 1 +4 by Rea et al. (2005). A possible geometry is a thick- 
ened or perhaps warped disc. Such a geometry in the case of Cyg 
X-3 has to be symmetrical to the compact object, not to cause any 
significant spectral changes related to the orbital phase, which are 
not observed. 

If the system inclination is lower, the observed strength of the 
Compton reflection component corresponds to a smaller value of R, 
and no corrections to the luminosities given in Table 2 are required. 
The upper limit on the inclination from the lack of eclipses for the 
two extreme binary models considered in SZ08 is 77° and 64° for a 
companion star mass of 5 and 50Mq, respectively. The higher value 
for the low mass of the companion star corresponds to a neutron star 
accretor, which we rule out. Since the star has an extreme stellar 



wind that is optically thick near the surface, the actual inclination 
has to be substantially smaller and is probably < 60°. 

Even if the presence of a strong reflection component can ex- 
plain the shape of the hard state, the unusually low electron temper- 
ature in the hard state, however, remains to be explained. The uni- 
versal cut-off at ~ 100-200 keV observed in both Galactic black 
holes (Zdziarski & Gierlinski 2004) and in Seyfert I type AGN 
(Zdziarski et al. 1996a), strongly suggests that the hard state is usu- 
ally characterised by electron temperatures of the order of 50-100 
keV. A tentative explanation is that the low electron temperature 
as well as the unusually high fraction of non-thermal electrons in 
the hard state is somehow a result of the turbulent environment with 
the accretion flow enshrouded in the strong Wolf-Rayet wind. If the 
accreted energy is transferred to the electrons in the form of accel- 
eration rather than thermal heating, this could be a sign of constant 
shocks and collisions between the wind and the accretion flow, as 
discussed in Hj08. A highly efficient acceleration process may also 
be responsible for the unusually strong radio emission in Cyg X-3, 
even if this can also be explained simply by the high density sur- 
rounding medium increasing the radiative efficiency of the jet. 

An alternative scenario that would result in a similar shape of 
the hard state spectrum is down- scattering of initially much higher 
temperature photons in a dense wind. The optical depth of such a 
wind has to be around t ~ 3 - 5 in order to down-scatter photons 
with initial energies of 50-100 keV down to the observed cut-off 
at ~ 20 keV. This high optical depth can probably not be achieved 
in the surrounding Wolf-Rayet wind, which was found by SZ08 to 
have an optical depth of ^ 1 . A dense disc wind instead of a thick- 
ened disc could however cause down- scattering instead of single re- 
flection in a similar geometry. While such a scenario could explain 
the low electron temperature, it would imply a much stronger ini- 
tial high-energy luminosity than the one calculated from our mod- 
els which would point to an even higher mass of the compact object 
than the one suggested here. 

6.4 The mass of the compact object 

Our results favour a massive black hole accretor in Cyg X-3 and 
the presence of a neutron star could be ruled out due to the high 
luminosity at the observed state transitions. A value of M = 30Mq 
corresponds to the state transition taking place at L = 0.03Le, in 
agreement with values observed in other sources and consistent 
with our present theoretical models for advective flows. We stress 
that our given values for the luminosities are the most conservative 
ones with respect to our models. Firstly, we do not correct for any 
net loss of photons by scattering in the wind due to orbital mod- 
ulation (Section 3.2). Secondly, we do not correct the luminosity 
for the part of the emission obscured from our line of sight by the 
reflector (Section 6.2). The combined effect of correcting for both 
these effects may increase the given luminosities by up to a factor 
of 1.7. 

A massive back hole is consistent with the results of Schmutz, 
Geballe & Schild (1996) who derive a mass for the compact object 
of 7-40Mq for a range of Wolf-Rayet masses from 5 to 20Mq and 
inclinations 30-90°, for the mass function of 2.3Mq. Limiting the 
inclination to < 60° and allowing a Wolf-Rayet mass between 5- 
50Mq, this mass function gives a mass of the compact object of 
9-60Mq. A black hole mass of 30Mq would imply a Wolf Rayet 
star mass of < 57Mq. A lower inclination allows for a lower Wolf- 
Rayet mass. 

Interestingly, recent results from mass measurements of the 
only two other Wolf-Rayet X-ray binaries discovered so far, IC 
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10 X-l (Prestwich et al. 2007, Silverman & Filippenko 2008) and 
NGC 300 X-l (Cai'pano et al. 2007a, b), suggest that the compact 
object in both cases are black holes and in at least one of them, IC 
10 X-l, the most probable mass is within 23-34Mq. 



7 CONCLUSIONS 

We have modelled Cyg X-3 in all spectral states as observed by 
RXTE, and for the first time derived unabsorbed spectral shapes 
and luminosities. Our modelling suggest that in 4 of the 5 spec- 
tral states, Compton reflection exceeds R = 2, corresponding to the 
maximum solid angle in an unobscured geometry, indicating that 
the emission from the innermost accretion flow is partly obscured 
(for i > 60°). We suggest that the peculiar shape of the spectrum 
in the hardest state of Cyg X-3 is caused by this strong Compton 
reflection on top of an intrinsic moderately hard state with a trun- 
cated disc, where the plasma electron temperature is much lower 
and the non-thermal electron fraction much higher than what is usu- 
ally observed in other systems. The low electron temperature could 
be due to down-scattering in a dense disc wind. The soft states are, 
corrected for absorption, very similar to those displayed by other 
sources like GRO J1655-40 and GRS 1915+105. 

We discuss the spectral evolution of Cyg X-3 and compare 
with other sources, using colour-colour and colour-luminosity di- 
agrams. In our interpretation of the intrinsic absorption corrected 
spectra, Cyg X-3 spans the whole range of spectral states as seen 
in other black hole sources. This is also consistent with its radio 
behaviour (Szostek et al. 2008). We find that the track traced out 
by Cyg X-3 in the colour-colour diagram is very different from that 
of the neutron star atoll sources but similar to that of black hole 
binaries, even if the Cyg X-3 colours do not enter the 'black hole 
only' regime as defined by DG03. Based on our most conservative 
estimates of the intrinsic luminosities, we study the spectral evo- 
lution as a function of Le, assuming three different masses for the 
compact object. We find that a neutron star accretor is not a likely 
candidate since not only the brightest states but also the transition 
to the hard state would then take place at luminosities close to Ed- 
dington, which is neither observed in other sources nor supported 
by our present theoretical understanding. For the state transition to 
take place at the expected L at 0.03Le, the compact object must be 
a black hole of M at 30M©, similar to the Wolf-Rayet black hole 
binaries recently discovered in other galaxies. 
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